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Abstract 

The evolution of two, non- inter acting two- level atoms immersed in a weakly coupled 
bath can be described by a refined, time coarse grained Markovian evolution, still 
preserving complete positivity. We find that this improved reduced dynamics is able 
to entangle the two atoms even when their internal frequencies are unequal, an effect 
which appears impossible in the standard weak coupling limit approach. We study in 
detail this phenomenon for an environment made of quantum fields. 



1 Introduction 

Independent, non- interacting atoms immersed in a common bath, represent an example of 
an open quantum system, i.e. a subsystem in interaction with a large environment p-J-[6j. 
Their time evolution can be given as a reduced dynamics, obtained by eliminating the degrees 
of freedom of the environment and by subsequently performing a Markovian, memoryless 
approximation, justified by the very rapid decay of correlations in the bath. The resulting 
atom evolution is irreversible and incorporates dissipative and noisy effects induced by the 
environment. 

In many physical instancies, the atoms can be treated in a non-relativistic approxima- 
tion, as independent finite-level systems, with negligibly small size. On the other hand, the 
environment can be modeled by a set of weakly coupled quantum fields, typically the elec- 
tromagnetic field, in a given temperature state, interacting with the atoms through a dipole 
type coupling [7j. Although ignoring the internal atom dynamics and the full vectorial struc- 
ture of the electromagnetic field, this simplified setting is nevertheless perfectly adequate for 
studying the behaviour of physical systems like ions in traps, atoms in optical cavities and 
fibers, impurities in phonon fields [Z]-[9]. 

The derivation of an acceptable subdynamics for the atoms is notoriously tricky [10], PQ-[3] 
and time evolutions that are not even positive have been adopted in the literature in order to 
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describe their physical properties [5], [TT] [12] . Instead, a physically consistent time evolution 
for the atom subsystem can be obtained through a suitable coarse grained procedure, within 
the weak coupling approximation [H [3], [13[ [HI EE] : the resulting subdynamics is described by 
a one parameter (= time) family of completely positive maps that form a quantum dynamical 
semigroup. 

We shall explicitly discuss below such a derivation for a subsystem composed by two, 
unequal, mutually non-interacting atoms. For simplicity we shall restrict the attention to 
two- level atoms weakly coupled to a collection of independent, free, massless scalar fields in 
3 + 1 space-time dimensions, assumed to be in a state at temperature T = 1/(5. 

As well known, the interaction with an environment usually leads to decoherence and 
noise, typical mixing enhancing phenomena. Therefore, one generally expects that quan- 
tum correlations that might have been present at the beginning between the two atoms be 
destroyed when they are immersed in the bath. 

However, an external environment can also mediate indirect interactions between otherwise 
totally decoupled subsystems and therefore a mean to correlate them. This phenomenon is 
generic on a short, microscopic time-scale where the dynamics is unitary and reversible; 
however, its persistence is not expected in general on longer time-scales where irreversible, 
dissipative and decohering effects described by Markovian master equations appear. 

Nevertheless, there are instances where purely dissipative, non-Hamiltonian contributions 
to the master equation can lead to entanglement generation [H)]-[23],[E]- This is due to 
the particular form of the (Kraus) operators appearing there, which couple indirectly {i.e. 
not dynamically) the two subsystems. This phenomenon has been established in the case 
of subsystems formed by two identical two- level systems [23]- [27] or harmonic oscillators 
[2"5| |2"§] evolving with a reduced dynamics obtained via the so-called weak coupling limit [30J . 
This technique is applicable when the time-scale over which the dissipative effects become 
visible is so large that the free dynamics of the subsystems can be effectively averaged out, 
thus eliminating very rapid oscillations; this typically occurs for environments with very fast 
decaying correlations [3], [5] . Nevertheless, it turns out that, in the case of atoms with 
unequal frequencies, this procedure of averaging out fast oscillations prevents generation of 
entanglement; in other terms, the environments for which the weak coupling limit procedure 
is justified are unable to correlate two atoms with different internal frequencies, while they 
do when the atoms are identical. 

In the following, we shall study in detail the conditions that allow the two otherwise inde- 
pendent, unequal atoms to become initially entangled through the action of the environment, 
when the weak coupling limit procedure is not applicableo Following [15], we shall instead 
derive a time evolution for the two atoms that allows a finite coarse-graining time interval: 
remarkably, it results still expressible in terms of a completely positive quantum dynamical 
semigroup, which reduces to the standard one obtained through the weak coupling limit in 
the limit of an infinitely large coarse-graining time interval. We shall see that in this refined 

1 A preliminary investigation on these topics in the case of two qubits weakly coupled to a Ohmic bath 
has been reported in [31]. 
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framework the sharp dependence on the atom frequencies of the entanglement capability 
of the environment mentioned above looks more like a mathematical artifact than a real 
physical effect. 

2 Two-atom reduced dynamics 

As explained above, we shall study the behaviour of a system composed by two, unequal two- 
level atoms, that start interacting at time t — with an environment made of a collection 
of independent, massless, scalar quantum fields at temperature We are interested 

in the evolution of the atoms as open quantum systems and not in the details of their 
internal dynamics; therefore, we shall model them, in a nonrelativistic way, as simple qubits, 
described in terms of a two-dimensional Hilbert space. 

In absence of any interaction with the external fields, the single atom internal dynamics 
can thus be taken to be driven by a generic 2x2 hamiltonian matrix. As a result, the total 
atom Hamiltonian Hs can be expressed as: 

3 

U _ pH 1 ) I rr(2) rr(«) _ U a> _ (a) _ , 9 m 

h s -h s +h s , H s - —n- a" = — 2_^n i a i , a-1,2, (1J 

i=i 

where a\ = cr, <8> 1 and op = 1 <g> t7j are the basis operators pertaining to the two different 
atoms, with Gi, i = 1,2,3 the Pauli matrices, rii the components of a unit vector, while uo a 
represent the gaps between the two energy eigenvalues of the two atoms. 

As mentioned in the introductory remarks, the coupling of the atoms with the external 
fields is assumed to be weak, so that the dipole approximation results appropriate [7] . In our 
simplifying settings that ignore spinorial indices, the interaction term can then be described 
by an hamiltonian H' that is linear in both atom and field variables: 

H' = j^{°? ®Uf {l) \+°? ) ®Uf {2) \) ■ (2) 
t=i 

The operators x) represent the set of external quantum fields, taken to be spinless 
and massless for simplicity. They evolve in time as free relativistic fields with a standard 
Hamiltonian Hq, [32] • The atoms are assumed to have a spatial extension described by the 
two functions f^ a \x), a = 1, 2. To be more specific, we shall choose for the atoms a common 
profile f(x) of spherically symmetric shape, with size e: 

m = _l m 2 , (3) 

* [|2| 2 + (£/2) 2 ] 

2 As we shall see in the following, the choice of modelling the environment in terms of relativistic quan- 
tum fields (at finite temperature) allows an analytic treatment of the reduced two-atom dynamics, without 
additional approximations besides the Born and Markov ones (see below). 
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and position the first atom at the origin of the reference frame, so that f^\x) = f(x), while 
the second is displaced by an amount £ with respect to it, f^ 2 '(x) = f(x+£). Since the atom- 
field interaction takes place on the whole region occupied by the atoms, the field operators 
entering the interaction Hamiltonian above are smeared over the atom spatial extension: 

Uf {a) } = jd 3 xf^(x)^{0,x) , a = l,2. (4) 

The total Hamiltonian H describing the complete system, the two atoms together with 
the external fields $j, can thus be written as 

H = Hs + Hq + XH' = H + XH' , (5) 

with A a small coupling constant H Through the standard Liouville-von Neumann equation, 
d t p tot (t) = —i[H, p tot (t)], it generates the evolution in time of the state of the total system, 
described in general by a density matrix p tot , starting at t — from the initial configuration: 
Ptot(O). 

We shall assume the atom and the fields to be initially prepared in an uncorrelated state, 
with the fields in the temperature state pp and the atoms in a generic initial state p(0), so 
that ptot(O) = p(0) ® pp. The reduced time evolution of the two atoms is then obtained 
by integrating over the unobserved field degrees of freedom and is formally given by the 
transformation map: p(0) t— > p(t) = Tr$[p to t(*)]- This map is in general very complicated, 
because of nonlinearities and memory effects; nevertheless, it can be approximated by a lin- 
ear, memoryless map when the coupling with the environment is small and its own internal 
dynamics is sufficiently fast Indeed, in such cases the details of the internal environ- 

ment dynamics result irrelevant, being the time scale of the subsystem evolution typically 
very long compared with the decay time of the correlations in the bath0 

In order to derive the equation obeyed by the reduced density matrix p(t) in the case at 
hand, it is convenient to work in the interaction representation 

p to t (t) = j tHo Ptot if) e- itH ° , (6) 

so that 



<9ptot(*) -x 



H'(t), p tot (t) , H'(t) = j m ° H> e- itH ° . (7) 



3 For simplicity, we have assumed that the environment couples with the same strength to the two atoms: 
this makes the following analytic derivation of the master equation more transparent, without compromising 
its generality. 

4 A discussion on the validity of this so-called Markovian approximation is reported in 13J . There, a 
non-Markovian weak coupling approximation of the reduced dynamics is also introduced; it leads to a two- 
parameter family of dynamical maps, with a time-dependent generator We stress that this approach is 
completely different from the one discussed below, which instead describes the reduced two-atom dynamics 
in terms of a Markovian, one parameter semigroup. In particular, while in [13] the standard weak coupling 
limit can be reached only in the asymptotic, long-time regime, in the treatment presented below it can always 
be obtained for any time by letting the coarse-graining parameter become large. 
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One then focuses on the changes of the reduced state p(t) = Tr$[p tot (t)] over a time interval 
At; by taking the trace over the field variables of the integrated version of the equation (JZ|) 
one gets (to lowest order in A): 

pa + At) - pa) i r t+At ds dp(s) 



At At J t ds 

2 rt+At rti 



— J dh J dt 2 Tr^H'ih), [H'(t 2 ), (*)]]) +0(A 4 ) .(8) 



One notices that the variation of p(t) starts to become relevant at order A 2 , i.e. on time 
scales of order r = X 2 t. Then, one can equivalently write: 

p(t + At) -p(t) 1 [ T+x2At ds dpjs/X*) ^ 



At At J T ds 

so that in the limit of small A (and finite At) one can readily approximate the r.h.s. of ([9]) 
with dtp{t). At this point, one further observes that the environment, containing an infinite 
number of degrees of freedom, is much larger than the subsytem immersed in it, so that its 
dynamics is hardly affected by its presence. It is therefore justified to replace in the double 
integral of (jHJ) the evolved total state ptot(t) with the product state p(t) £g> pp, taking the 
initial state pp as a reference state for the bath [I]-[3j. 

Returning to the Schrodinger representation, one finally gets the following linear, Marko- 
vian master equation for the two-atom state pit): 

^- = -i[H s ,p]+V[pit)], (10) 

where the bath-dependent contribution £>[p(t)] contains both an Hamiltonian and a dissipa- 
tive term 

V[p(t)} = -i[H 12 , pit)] + jC[p(t)} , (11) 

with 

?a 2 r At r At / \ 

#12 = ^ t J o ds ^ J Q ds * 6 ( s i ~ s 2) Tr $ (p^ [H'( Sl ), H'(s 2 )]) , (12) 



L(p{t)®pp)L--{L 2 ,pit)®p p } 



At 



dsH'is) , (13) 



o 



the curly brackets representing the anticommutator, while 6*(s) is the step function. 

It is important to observe that, for any interval At, the master equation (|T0|) - (|T3|) generates 
a quantum dynamical semigroup of completely positive maps. Indeed, the generator in the 
r.h.s. of fllip . besides the Hamiltonian piece, contains a dissipative term which turns out to be 
itself completely positive, being the composition of two completely positive maps, the trace 
over the environment degrees of freedom and a linear operator on the total system, written in 
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canonical Stinespring form [331 EHE]- Notice that, on the contrary, in the usual weak coupling 
limit approach to the derivation of a Markovian master equation, complete positivity is 
ensured by an ergodic average prescription, that, as mentioned in the Introduction, eliminates 
fast oscillating terms [TQ1 ED]; in the present formalism, this corresponds to letting the time 
coarse-graining parameter At going to infinity: it is thus applicable only to environments 
with sharp decaying correlations. In the following, we shall instead keep At finite and 
consider therefore more general situations. 



3 Master equation 

For the case at hand, a more explicit expression for the generator in fill I) can be obtained 
by recalling (J2J) and (HI) . Indeed, after straightforward manipulations, the master equation 
driving the dissipative dynamics of the two atoms state takes the following Kossakowski- 
Lindblad form [351 [36] 

iW) --i[H eS: p{t)]+C[p{t)\ , (14) 



with 



and 



dt 



EE". 

a,/3=l i ,j=l 



(a/3) (a) (/3) 
ij i 3 



m = E E c l 



(a/3) 



a,/3=l i,j=l 



03) H 



(a) iP) 



(15) 



(16) 



The coefficients of the Kossakowski matrix C\j and of the effective Hamiltonian H e g are 
determined by the field correlation functions in the thermal state pp. 

G\f\t -t') = J d*xd 5 yf^\x) f^(y) (^(t,x)^(t',y)) , (17) 



through their Fourier, 



and Hilbert transform, 



dte lzt G[f\t) 



dtsign(t)e izt Glf \t) = - I dw ij ' ' ' 



w — z 



(19) 



respectively (P indicates principle value). 
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More specifically, one finds that the Kossakowski matrix reads: 



3 

c {a/3) = x 2 V V e i(e<»*+e'up)&t/2 ^,(0 jjti') 

ij / j / j i hi i Ij 

£,£'={+ -,0} k,l=l 

x At r H, , U A Sin ^ - ^ At / 2 } sin [(go; + fag) At/2] 

^•f = rij % , = - (% - rij ± ie ijk n k ) , (21) 

are the components of auxiliary three-dimensional tensors, giving the free evolution of the 
atom operators: 

3 



where 



' (t) = e m *a\ a) e~ ltHs = ^ °\ 

?={+,- 0} j=l 



The 6x6 matrix C|- turns out to be non-negative, since, as already mentioned, the 
evolution generated by Q) is completely positivejj An expression similar to the one in (1201) 
holds also for in (IT51) . with G^ ( u ) replaced by K k ^\ui). 

For simplicity, the fields giving rise to the environment are taken to be independent and 
further assumed to obey a free evolution; in this case, one finds: 

(Mx)$M) = ^[M^MPp] = ^G{x - y) , (22) 

where G(x — y) is the standard four- dimensional Wightmann function for a single relativistic 
scalar field in a state at inverse temperature (3 [32], that, with the usual ie prescription, can 
be written as: 

G(x) = J -0^ 9(k°) 5(k 2 ) [(1 + Af(k )) e~ ik - x + Af(k°) e ik ^ e~ £k ° , (23) 

where 

Although the ie prescription, assuring the convergence of the integral in f[2"3l . originates 
from causality requirements, in the present setting it can be related to the finite size of the 
two atoms. Indeed, the correlations in (TTTjl actually involve the Fourier transform f(k) = 
J d 3 x e lk ' x f(x) of the shape function f(x) in (j3J); it can be easily computed to be f(k) = 
e -|fc|e/2_ inserting it back in (|17|) . this contribution can be conveniently attached to the 



5 On the other hand, let us remark that direct use of the standard second order perturbative approximation 
(e.g. see [E1IH]) often leads to physically inconsistent results [S1[TTJ[T2], gi vm § a finite time evolution for p(t) 
that in general does not preserve the positivity of probabilities. 
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definition of the Wightmann function G(x), so that the integrand in (123]) gets an extra e £k ° 
overall factor. 

Using (1251) and ([21]), the Fourier transform in (jTSJ) can now be explicitly evaluated; taking 
for simplicity the limit of pointlike atoms, (the size e can be taken to vanish since it does 
not play any more the role of a regularization parameter), one gets: 



"I J 



LO 



8 l3 



with: 



g^(co) = g^\cu) 



1 



10) 



uo 



(25) 



2tt 1 - e-^ ' 
1 uo sm(£uo) 



2tt 1 - e-P" £uo 



(26) 



where £ denotes the modulus of the displacement vector £; then, recalling (fT9~j) . for the Hilbert 
transform one similarly finds: 



7TZ 



G?U> 



5^1 



u> — 2; 



(27) 



With these results and taking into account that Ylik^H ^Ij = "(£ + the Kos- 
sakowski matrix takes the more explicit form: 



fc=0 



where 



(a/3) 



with 



r(«/3) 



while 



At 

Air 



duo 



g (al3 \uo)±g (al3 \-uj) 



'0 



sin [(uo - u a )At/2] sin [(uo - u p )At/2] 
(uo - uo a )At/2 (uo - up) At/2 

sin {uoAt/2) 



g^)(uo) + g^\-uo) 



1 2 



uoAt/2 



(28) 



(29) 



(30) 



(31) 



Only the following combinations g±^\uo) = (?( a/3 )(u;) ± £( a/3 )( — cj) actually occur in the 
previous integrals, and from the explicit expressions in (12"6"1) one obtains: 



,(11) 



,(11) 



r {22) _ 00 
Z7T 



1 + e 



-0u> 



1 - er&> 



,(22) 



2r7 



,(12) 



,(12) 



y+ ~2n 



1 + e 



1 - e-^ 
(21) _ uo sin(uo£) 
" " 2^ oj£ ' 



sin(w 



(32) 
(33) 



they contain the dependence on the bath temperature 1//3 and on the separation £ between 
the two atoms. Because of the presence of the Boltzmann factors, the integrals in 



( |30|) . (13T]) can not in general be expressed in terms of elementary functions. However, in the 
case of a bath at high temperature (i.e. for small (3), the square bracket in (132]) behaves as 
2/ (3u and the above integrals can be explicitly evaluated (see the Appendix). In the physical 
situation for which £ < AtJ^I one finds: 



r(°?0) 



r(«/3) 



sin (£ua/2) 



nf3uj a p At I £u a /2 



sm 



sin Uup/2) 
H ; — — - sm 



eup/2 



(u a (l-£/At) — 



1 



Tri UJ a fi At 

1 



sm 



u a/3 (At - 



sm 



£(u a + up) 



4tt(3 \ 2 At 



(34) 
(35) 
(36) 



Inserting these results back in (129|) and (l3Ti) , one finally obtains the explicit expression for 
the Kossakowski matrix in the large temperature limit 

Coming now to the Hamiltonian contribution to the master equation, one sees that the 
effective Hamiltonian H e g in ([151) can be split into two parts, H eS = H s + the first is 

just a renormalization of the starting system Hamiltonian, while the second one represents 
an environment induced direct coupling term for the two atoms. The term Hs has the same 
form as the Hamiltonian in ([!]) but with redefined frequencies 



LO n — l- 



At 

'27T 



UJ 



sin Uuj — u a ) At/2] 
(u - u a )At/2 



(37) 



Recalling the definition of K,^ aa \uj) in (1271) . one sees that it can be split as: 



JC (aa) { 



2lX 2 l 



P / dz 



z — UJ 



z 



1 



1 



+P dz , , , 

/ 1 — & 3Z \ z + UJ z — UJ 



(38) 



6 This condition assures that the two atoms actually feel the presence of the quantum fields; indeed, due 
to relativistic causality |32j . the fields would not be able to interact with the atoms in the time interval At 
if they were to far apart. 

7 Since this is an approximated result, positivity of the matrix is not a priori guaranteed and should be 
formally imposed in order to preserve the properties of the exact expression (1201) . In particular, positivity of 
the two diagonal submatrices c\" a \ requires (3oj a /2 < 1, which are satisfied by the requirement of /3 small. 
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into a vacuum and a temperature-dependent piece. Although not expressible in terms of 
simple functions, the temperature dependent second term is a finite, odd function of uj\ on 
the contrary, the remaining, vacuum contribution in fl38l) results divergent, and therefore 
so are the shifted frequencies uj a . As a consequence, the definition of effective Hamiltonian 
H e g requires the introduction of a suitable cutoff and a renormalization procedure. This is 
not a surprise: the appearance of the divergences is due to the non-relativistic treatment 
of the two-level atoms, while any sensible calculation of energy shifts would have required 
the use of quantum field theory techniques [7j. In order to make H e R well defined we follow 
a simple prescription: perform a suitable temperature independent subtraction, so that the 
expressions in (T5TT) reproduce the correct quantum field theory result, obtained by considering 
the external fields in the vacuum state. 

( 12) 

The induced two-atom interaction term if eff can instead be expressed as 



where 



(12) 



H 



with 



Jo 



COS 



H, 



(12) 



cff 



E" 

*>J=1 



toi 2 At 



(k+sin 



2_^e ijk n k \J + +\Jq-cos — - — J+ I 
- fc=o ' ^ L J / 



-i—— I doj 

47T 



oo 
oo 



-i— - I (he 

87T 



to 



wn[(u-ui)At/2] sin [(u - u 2 ) At/2] 

(uj - uj^At/2 (uj - uJ 2 )At/2 

sin (uAt/2) 12 
uAt/2 



(39) 



riiUj , (40) 



(41) 
(42) 



Also K,^ {uj) can be split as in (13"5|) into a temperature dependent term, odd in u>, and 
a vacuum piece. Clearly, only this second contribution enters the above integrals J + ; it 
is finite (for non vanishing atom separation) and with the help of fl32|) can be explicitly 
computed: 



Inserting this result in (|4"T|) and 
1 



P 
2~kH 



dz 



sin £z 



1 cos £uj 



(43) 



z + uj iz 2ni £ 
, one finally obtains, again for I < At (see Appendix): 



J4 



2n£oj 12 At 



cos 



~(coi + uj 2 )£' 


sin 


"a;i 2 (At-^) 


to 




2 



Jn 



If 1 _I 



(44) 



We are now ready to discuss the entanglement properties of the time evolution generated by 
the master equation fn4l) - ffT6l) . 
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4 Environment entanglement generation 



In order to study the entanglement power of the thermal bath made of free quantum fields, 
we shall focus on the small t behaviour of the dynamics generated by the equation (JT3J); our 
aim is to investigate whether the two independent atoms can get entangled by the action of 
the environment in which they are immersed at the beginning of their dissipative evolution. 
Then, without loss of generality, one can limit the considerations to pure, separable initial 
states, and therefore take: 

p(P) = \<p)(<p\®\i{>)(i!>\ , (45) 

with given single atom states; indeed, if the environment is unable to create entan- 

glement out of pure states, it will certainly not correlate their mixtures. 

Since we are dealing with a couple of two-level systems, one can use partial transposition 
as a criterion for entanglement creation [371 138] . More precisely, the environment is able to 
create quantum correlations between the two atoms if and only if the operation of partial 
transposition spoils the positivity of the state pit). 

The presence of negative eigenvalues in the partially transposed reduced density matrix 
pit) can be ascertained by looking at the sign of the average 

A(t) = ( X \ P(t) \X) , (46) 

with \x) a four-dimensional vector. Indeed, choose \x) to be orthogonal to \ip) ® \ip), so that 
the above average initially banishes, -4(0) = 0. Then, the two atoms, initially prepared in 
a state p(0) = p(0) as in fj4"5l) . will surely become entangled if \x) can be further chosen so 
that dtA{0) < 0. From this condition, a simple test for entanglement creation involving the 
elements of the Kossakowski matrix (1281) and of the effective interaction Hamiltonian (1391) 
can then be extracted [2UES]. It explicitly reads: 

(u\C^\u) (v\{C^) T \v) < \(u\Ke{C^+iH^ ] )\v)\ 2 , (47) 

where T means matrix transposition; the three-dimensional vectors \u) and \v) contain the 
information about the starting factorized state (I45I) : in fact, their components can be ex- 
pressed as: 

Ui = ((p\<7i\ip ± ) , Vi = (ipx\ai\ip) , (48) 

where \<p±) and \ip±) are the orthonormal complement to the initial atom states \(p) and 
respectively. Therefore, the external quantum fields will be able to initially entangle the 
two atoms evolving with the Markovian dynamics generated by (fT4"|) if there exists an initial 
state of the form fT4oT) for which the inequality (I4"7j) is satisfied. 

In order to obtain a manageable expression for it, we first note that, without loss of 
generality, the unit vector n that defines the atom Hamiltonian in (pQ) can be oriented along 
the third axis. Further, as initial atom state we shall choose p(0) = |— )(— | ® |+)(+|, 
constructed out of the eigenstates of the single atom Hamiltonian, er 3 |±) = ±|±). As a 
consequence, recalling (JHJ), one finds that the three-dimensional vector \u) has components 
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Ui = {l,—«,0}, and further t>j = U{. Then, using the explicit expressions for the elements 
of the Kossakowski matrix C^- and of the induced interaction Hamiltonian the 
inequality (1471) reduces to 

~ ^t) i 1 + ^t) < ^[v+y +4 ( j + )i • < 49 > 

Notice that the l.h.s. of this expression is positive, since as discussed in the previous Section, 
complete positivity requires j3u a /2 < 1. 

As remarked at the end of Section 2, the parameter At identifies the time scale over which 
the presence of the environment is felt by the system of the two atoms; clearly the weaker the 
coupling with the environment is, the longer one needs to wait for the bath induced effects 
to become apparent. 

Let us first discuss the standard weak coupling limit approximation; in this case, one 
actually let the coupling constant A to approach zero, so that changes in the two-atom density 
matrix become visible only for infinitely large At. In this limit however, the two integrals on 
the r.h.s. of (B9i) become vanishingly small Jf| Thus, for atoms with unequal frequencies, the 
inequality can never be satisfied, and thus no entanglement is generated. On the contrary, 
when the two frequencies coincide, u>i = u 2 = u, the condition (j4"9~j) becomes: 



2 



sin [u 



2 



4 



COS [LO 



(50) 



This result generalizes the one discussed in [26] , where the contribution of the environment 
induced interaction Hamiltonian (the second term in the r.h.s. of fj49l) ) was neglected. In 
particular, one sees that, in this case, for any given (small) inverse temperature /3, there is 
always an atom separation i below which the inequality (15"0"|) is satisfied, and therefore en- 
tanglement created between the two atoms. This phenomenon is forbidden only for infinitely 
large separation or infinitely large temperature, in which case the environment induced de- 
coherence and noisy effects dominate. 

The sharp dependence of the entanglement capability of the environment on the atom 
frequencies in the weak coupling limit approach is however striking; it originates in the 
elimination of fast oscillating terms in the reduced two-atom dynamics through an ergodic 
average, a procedure that is justified only in the limit of a vanishing A and very fast decay 
correlations in the environment. 

Instead, if the coupling of the atoms to the bath is weak, but not infinitesimally small, 
environment induced changes in the atom density matrix pit) can be seen on finite time 
intervals At. In this case, it is the full condition (|49i) that regulates the entanglement 
capability of the thermal bath. One can check that indeed this inequality can be satisfied 
even for u\ ^ uj 2 , and therefore that a bath made of thermal quantum fields can correlate 
two unequal atoms. 



8 Indeed, in the limit At — > oo both integrals vanish since, for u>i ^ ui2 , the two functions 
sin[(w — uji)At/2]/[(ui — uii)At/2] and sin[(w — w 2 ) At/2]/[(w — ui 2 )At/2] have disjoint supports. 
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In order to show this, let us first note that the Hamiltonian contribution in (j49[) . being 
positive, can only enhance entanglement generation; this is the result of the hermiticity of 
the induced coupling term in ([39]) rl One can therefore limit the considerations to 

a simpler inequality, in which the term (J+) 2 is neglected; when this reduced condition is 
satisfied, also the full one in will clearly be. Recalling (1511) . and keeping for simplicity 
only first order terms in £, the condition for environment assisted entanglement generation 
reduces to: 



1-^1 (l + ^l\< 



sin (u l2 At/2) 
U12 At/2 



sin (u 12 At) 
At J u 12 At ' 1 ' 



For given bath temperature and atom frequencies, this condition is satisfied provided a 
sufficiently small time interval At results from the coupling with the environment. Further, 
the smaller the atom separation is, the easier the condition of (I5T|) will be met. 

Finally, note that, in contrast to situation encountered in the weak coupling limit approx- 
imation, here there is no sharp change between the regime of entanglement generation and 
the one of solely decoherence; the transition is smoothly regulated by the coarse graining 
parameter At, i.e. ultimately by the strength of the coupling of the atoms to the environ- 
ment. 



5 Discussion 

We have seen that two atoms, prepared initially in a separable state, can get entangled as a 
result of their independent interaction with a common bath made of thermal quantum fields 
even when their internal frequencies are unequal. This result is based on a novel Markovian 
approximation of the reduced atom dynamics, that allows an explicit dependence on the time 
scale At, measuring the interval over which the atoms feel the presence of the environment. 

This conclusion contrasts with the one obtained through the usual weak coupling limit 
approach to the atom reduced dynamics; in that case, the entanglement power of the external 
environment is reduced to zero for atoms with unequal frequencies as a consequence of the 
procedure of taking the ergodic average. In the light of the results presented in the previous 
Section, this conclusion appears however a mathematical artifact, originating in letting A 
to zero and At to infinity] 10 ! conditions hardly met in actual physical situations. Instead, 
for weakly coupled baths with finite At, environment assisted entanglement generation is 
always allowed, and can be controlled through the external parameters, the bath inverse 
temperature (3 and the atom spatial separation £. 

In the high temperature case (i.e. (3 small) and arbitrary At, we have explicitly shown 
that this conclusion holds because of the condition fl5T|) . Similarly, in situations allowing a 

9 In general, the dissipative and Hamiltonian contributions in the r.h.s. of (|47[) can destructively interfere 
making the inequality harder to be satisfied and thus reducing the entanglement power of the environment. 

10 Because of the Riemann-Lebegue lemma, the ergodic average, on which the so-called "rotating wave 
approximation" is based, is strictly justified only in the limit A — > and At — > oo. 
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large, but finite At, a different approximation of the full entanglement condition (j4Tj) can be 
given; it can be obtained using techniques and procedures analogous to the ones discussed in 
the previous Sections. Neglecting again the Hamiltonian contribution, one explicitly finds: 



l-i?i 



where 



Ra 



1 


"sin (u 12 At/2)' 


2 


'fun R 2 \ 


4 


lu 12 At/2 




\002 Rl) 



1/2 



1 + e-P^ 



sm 



a: 



UJ 2 Rl 

uji R 2 



1,2 



1/2 



So 



(52) 



(53) 



It is a further generalization of the condition discussed in [26J in the case of identical atoms, 
to which it reduces for uoi = uj 2 and At infinite. Although valid only for large (but finite) 
At, it can always be satisfied with suitably chosen f3 and £. In particular, (1521) is always 
true in the zero temperature case, i.e. in the limit (3 — > oo; in other words, a bath made of 
quantum fields in the vacuum state is always able to generate entanglement, for any finite 
spatial separation of the two atoms. 

All the above considerations are based on the condition fl4"Tj) for entanglement enhance- 
ment; when satisfied, it assures that quantum correlations among the two atoms are gen- 
erated as soon as t > 0. It is however unable to determine the fate of this quantum cor- 
relations as time increases and in particular in the asymptotically long time regime. On 
general grounds, one expects that the effects of decoherence and dissipation that counteract 
entanglement production be dominant at large times, so that no entanglement is left in the 
end. There are however instances in which the entanglement generated at the beginning of 
the evolution persists also for asymptotically long times [271 [28| [5] . In order to fully clarify 
this situation, a complete study and classification of the set of the equilibrium states of the 
refined master equation f|TU]) - f|T3^ is necessary!"] Work on this topic is presently in progress 
and will be reported elsewhere. 



Appendix 

We indicate here how to compute the integrals that appear in the expressions of the Kos- 
sakowski matrix C\^\ (I3"U1) and (I3T|) . and in that of the effective Hamiltonian interaction 

term H eS , fj4"Tl) and fj4~2l) . In the high temperature case, the explicit evaluation of fl30l) 

11 Only partial results on the classification the equilibrium states of completely positive quantum dynamical 
semigroups have been so far obtained [H 0D] • 
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involves the computation of integrals of the following two types: 

sin(a; — a) sin(x — b) 



h 
h 



dx sinfcxl 



+00 



dx 



x — a x — b 
sin(cx) sin(x — a) sin(x — b) 



x 



x — a 



x — b 



(54) 
(55) 



with a, b, c positive constants. By decomposing the products of trigonometric functions in 
terms of linear combinations of sines and cosines, one can split e.g. I\ into the sum of three 
simpler integrals: 

h = I + 1(c) - J(-c) , (56) 



with 



cos(a — b) 



dx 



sin cx 



[x — a)(x — b) 



1(c) 



+00 



dx 



sin [(2 



c)x — a 



"] 



(x — a)(x — b) 



(57) 

By first changing the integration variable to y = (2 — c)x — a — b in 1(c), with c < 2, and 
then reducing the denominators in partial fractions in both integrands, one can express Iq 
and 1(c) as combinations of the following integral (e.g. see |41j): 



+00 



dx 



sm(ax) 
x + z 



n cos(az) , a > . 



(5? 



Explicitly, one finds: 



h = 
He) 



— 7T- 



cos(a — b) 



a 

7i 



4(a - b) 

so that, recalling (|56|) . one finally obtains: 



sin [(a + b)c/2] sin [(a - b)c/2] , 
^ cos [a(c — 1) + b\ — cos [a + b(c — 1)] 



(59) 
(60) 



Ii = 7rsin 



(a + b)c 



sin [(a-b)(l - c/2)] 
(a -b) 



(61) 



This result holds for c < 2; when c > 2, one is forced to use a different integration variable in 
the expression of 1(c) in (157|) . y' = (c — 2)x + a + b, and as a result ends up with a vanishing 
value for I\. As a function of the parameter c, the integral I\ is however continuous, since 
the expression in (I6T|) also vanishes at the boundary point c = 2. 
From the result flBTl). one further obtains: 



lim ' 

c^o V c 



sin(x — a) sin(x — b) 



71 



x — a 



x 



a + b\ sin(a — b) 



(62) 



15 



The integral I2 in fl55|) can be evaluated using similar manipulations. When c < 2, one 
explicitly finds 



7r / sin(ac/2) 



a — 



sin [a(l - c/2) - 6] + 



sin(6c/2) 



sin [a — 6(1 — c/2 



(63) 



while for c > 2, a simpler expression holds: 



h = 7T 



sin a sin b 



(64) 



Here again one sees that I2 is a continuous function of c, since the expression in (163]) reduces 
to the one in (1641) at the boundary value c = 2. Further, from the expression in (1631) . one 
easily obtains the following limiting results: 



lim I — 



c^O \ C 



and similarly, 



lim I 2 



sin(x — a) sin(x — b) sin(a — b) 

(IX : = 7T ; 

„ x — a x — b a — b 



-00 / \ / • \ 2 

Sill CI / SIM 



dx 



x 



X 



TTCll-- 



(65) 



(66) 



( 12) 

The integrals appearing in the evaluation of the Hamiltonian contribution H^ s can instead 
be all reduced to expressions of the form: 



J 



+ °° sin(a; — a) sin(x — b) 

dx cos(cx) 

^ x - a 



x — b 



(67) 



With the help of manipulations similar to the one used above, J can be reduced to combi- 
nations of the following integral [41] : 



+00 



dx 



cos[a x) 
x + z 



n sin(az) , a > 



(68) 



When c < 2, the integral in (]67l) can be cast in the following form: 



J = IT COS 



(a + b)c 



sin [(a-6)(l-c/2)] 



(a -b) 



(69) 



while it vanishes for c > 2. In the limit of vanishing a and b, it reduces to J — 7r(l — c/2). 
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